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q , Abstract. We present QCD prediction of double-spin asymmetries (Att) in transversely polarized 

Drell-Yan process at small transverse momentum Qt of dimuon. Resummation of large logarithmic 
^"^ | corrections, relevant in small Qj region, is performed up to next-to-leading logarithmic (NLL) 

accuracy. Att at RHIC, J-PARC and GSI are studied numerically in the corresponding kinematic 
regions. We show that the large Att is obtained for small Qt and moderate energies. 
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The transversity 8q(x) is one of the twist-2 distribution functions of nucleon, which 
represents the distribution of transversly polarized quarks inside transversly polarized 
\q . nucleon [1]. Although it is theoretically as important as other twist-2 distribution func- 

tions such as the unpolarized parton distributions (q(x),g(x)) and the helicity distribu- 
tions (Aq(x), Ag(x)), very little has been known of it so far. This is because 8q(x) is a 
chiral-odd function, and should be always accompanied with another chiral-odd function 
Qh| in physical observable, and therefore cannot be measured in inclusive DIS. Transversely 

polarized Drell-Yan (tDY) process, p^ p^ — > £ + £~X, p^ p^ — > £ + £~X, is one of the pro- 
cesses where we can access the transversity by measuring the double-spin asymmetry: 

Att = = ^TT+rfgU • Since tDY is inclusive in the final state, it is in principle the 
cleanest process to access 8q(x). However, at the RHIC-Spin experiments, the asymme- 
tries are likely to be quite small [0]. This comes from the fact that the DY process at 
pp collider probes the sea distributions which are likely to be small for the transversity. 
Moreover, the rapid growth of unpolarized sea distributions enhances the denominator 
of Att at low-x which is typically probed by tDY at RHIC. On the other hand, much 
larger asymmetries are expected at the proposed spin experiments at J-PARC and GSI, 
which are to be performed at lower energies |3|,|4|,l5|]. 

In this work, we study the asymmetries in Qt distribution of dimuon, especially at 
small Qt where the bulk of dimuon is produced. At small Qt, the cross section is not 
described correctly by the fixed order calculations since large logarithmic corrections 
such as Oc' s 1 log m (Q 2 /Q T )/Q T (m<2n— 1) appear at each order of perturbation series. 
These so-called "recoil logs" come from soft gluon emissions, and have to be resummed 
to all orders in a s to make a reliable prediction of the cross section at small Qt- 
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The resummation is carried out in the impact parameter b space, conjugate to Qt 
space, to enforce transverse-momentum conservation, and the resummed cross section 
is expressed as the Fourier transform back to the Qt space. Here we perform the 
resummation at the next-to leading logarithmic (NLL) accuracy, which corresponds to 
adding up the terms with m = In — 1 and In — 2, respectively, for all n. At this level, 
the "resummed part" of the spin-dependent cross section of tDY, differential in invariant 
mass Q, transverse momentum Qt and rapidity y of dimuon, and in the azimuthal angle 
<j) of one of the outgoing leptons is given by 10] (y/S is the CM energy of hadron system), 

A / Jg 7 NLL = cos(24>)— °^Yef f ' db-J (bQ T )e s ^ (1) 
dQ 2 dQ 2 dyd(j) V Y) 3N C SQ 2 ^ 'Jo 2 UV ^ J 

(C qq ® HO U, §) (Cqq ® Sqt) U, ^) + (X? - *g) 



Here Jo(bQr) is a Bessel function, bo = 2e r£ with Je the Euler constant. The large 
logarithmic corrections are resummed into the Sudakov factor e s ( b, @) with S(b, Q) = 

— f% 2 (dK 2 /k 2 ) {In A q (a s (K)) + B q (a s (ic))}. Perturbatively calculable quantities, 

i.e., A q , B q and the coefficient functions C qq (z) and C qq (z), are expressed as power series 
in a s , and the their explicit forms necessary at the NLL accuracy are found in Ref.|0]. 

2 = (Q/y/S)e ±y , and 8qi{x, }l 2 ) is the transversity of z-th flavour quark at the MS scale 
jU. The singularity in /^-integration, due to Landau pole in (X s (k), is taken care of by 
"contour deformation method" introduced in [8]. Correspondingly, the nonperturbative 

effects are included by the replacement e s ^ -> e s ^ Q ^ gNpb2 in © Q, H, @|, with 
a non-perturbative parameter g^p, which (roughly speaking) parameterizes the intrinsic 
/cr-distribution of quarks inside nucleon. Then we combine the resummed part (0Q), which 
embodies the logarithmically enhanced contributions for small Qt to all orders, with 
the "residual part" of the fixed-order cross section, which is not associated with such 
logarithmic enhancement. To perform this consitently, we need the leading order (LO) 
tDY cross section at finite Qt', this is of 0(a s ) and is obtained as QCD prediction at large 
Qt fzfa- The matching of the NLL formula (OQ) with the corresponding component in the 
LO cross section is performed at intermediate Qt following the formulation of Ref.|@], 
to ensure no double counting for all Qt, and we obtain the "NLL+LO" prediction of 
the tDY cross section, ATdo NhL+h0 /(dQ 2 dQjdyd^)), which has a uniform accuracy 
over the entire range of Qt QVD- The "NLL+LO" cross section of unpolarized Drell-Yan 
process, do NLL+h0 J(dQ 2 dQldyd^), is obtained in the same way, utilizing the results 
in the literature flSfl]. 

Taking the ratio of these cross sections, we obtain the double-spin asymmetries: 

A TT = [A r d(T NLL+LO /dQ 2 dQj-dyd(j)] /[da NLL+LO /dQ 2 dQ 2 T dyd§] . (2) 

In order to calculate Att numerically, we need to assume a model of the transversity 
8qi(x) for the numerator. Here we take a model used in [0], which saturates the Soffer 



bound 111 111 as dqi{x, jUq) = [qt(x, jUq) +Aqt(x, jUq)]/2 at the low input scale jUq ~ 0.6GeV 



and is evolved to the higher \i 2 with NLO DGLAP kernel [Q3]; as the inputs, we 
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FIGURE 1. A T t at NLL+LO accuracy in pp collision at RHIC kinematics: Vs = 200 GeV, Q 
3,5,8, 15 GeV and <j) = withy = (left panel) andy = 2 (right panel). 
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FIGURE 2. A TT at NLL+LO accuracy. Left: y/S = 10 GeV, Q = 2,3,4 GeV, = 0, and y = for 
pp collision at J-PARC kinematics. Right: ^/S = 14.5 GeV, Q = 2,3,4,6 GeV, <f> = 0, and y = for pp 
collision at GSI kinematics. 



use GRV98 [13] for /xj) and GRSV2000 [14] for A^(x,/^). Correspondingly, 
the GRV98 distributions are used for calculating the denominator of ©. The non- 
perturbative parameter g^p are taken to be common in the numerator and the denom- 
inator of ©, and we use g^p = 0.5 GeV 2 which is consistent with the result of lfl5ll .FI 

Fig. 1 shows A TT in pp collision at RHIC kinematics: ^/S = 200 GeV, 2 = 3,5, 8, 15 
GeV and 0=0, with y = (left panel) and y = 2 (right panel). In this case, A jj are 
about 4-8% and rather flat in the small Qj region where the resummed cross sections 
are dominant in both the numerator and denominator of ©. We have smaller Ajt for 
smaller Q due to the growth of the sea distributions in the denominator of © at small 
x. We obtain slightly larger Att for y — 2 compared with the y = case, and it appears 
that generically the v-dependence of Att is small. 

The left panel of Fig. 2 is same as Fig. 1, but at J-PARC kinematics: = 10 GeV, 



2 The ^/"-dependence of the polarized and unpolarized cross sections almost cancels in Att of <[2j> in the 
range g^/. = 0.3-0.8 GeV 2 . 



2 = 2,3,4 GeV, = and y = 0. In this kinematics, the parton distributions at medium 
x ( X 12 = 0-2-0.4) are probed, so that the transversity distributions in the numerator 
of d2l> are larger than the RHIC case while the growth of the sea distributions in the 
denominator is not significant. Therefore we obtain much larger Ajj ~ 15% than the 
RHIC case of Fig. I. Att are again flat as functions of Qt, and the dependence on Q is 
also weak. 

The right panel of Fig. 2 shows Ajt in pp collision at GSI kinematics: = 14.5 
GeV, Q = 2, 3, 4, 6 GeV, 0=0 and y = 0, where Att are dominated by valence distri- 
butions at medium x. The largest Att of 15-30% are obtained in this case. The results 
are extremely flat as functions of Qt- Integrating the numerator and the denominator of 
© over Q T , we reproduced the NLO asymmetries given by Barone et al.yfl. We also 
calculated the v-dependence of the results in Fig. 2, and it turns out to be small for both 
J-PARC and GSI cases. 

To summarize, we have calculated the double-spin asymmetries Att for small Qt 
DY pair production in high-energy pp collisions at RHIC and in moderate-energy pp 
and pp collisions at J-PARC and GSI. Our results demonstrate that Att reach a finite 
value at Qt = through the flat behavior in the small Qt region, reflecting that the 
soft gluon resummation to the NLL level has universal structure for the polarized and 
unpolarized DY [7], and also reveal that Att are large enough to be experimentally 
measured, especially in moderate energies. 
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